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Abstract 

The structure of heavy hyperhydrogen ^H is studied within the framework of a tnnA four-body 

cluster model. Interactions among the constituent subunits are determined so as to reproduce 

reasonably well the observed low-energy properties of the tn, tA and tnn subsystems. As long as 
O 
^ ' we reproduce the energy and width of ^H within the error bar, the ground state of ^H is obtained 



as a resonant state. 



I. INTRODUCTION 

One of the research goals in hypernuclear physics is to study new dynamical features 
by injecting a A particle into a nucleus. For example, since there is no Pauli principle 
between nucleons and a A particle, the A participation gives rise to more bound states and 
significant contraction of nuclear cores, especially in light systems. Such a dynamical change 
in light hypernuclei has been studied mostly in systems composed of a stable nucleus and an 



attached A particle 



QE 



2| . In light nuclei near the neutron drip line, there have been observed 
interesting phenomena concerning neutron halos. A corresponding subject in hypernuclear 
physics is to focus on structures of neutron-rich A hypernuclei. If a A particle is added 
to such nuclei with weakly-bound neutrons or resonant ones, a resultant hypernucleus will 
become more stable against neutron decay. Thanks to this gluelike role of an attached A 
particle, there is a new chance to produce a hypernuclear neutron- (proton-) state, if the 
core nucleus has an unbound (resonant) nucleon state with an appropriate energy above 
the nucleon-decay threshold. Another interest is to extract information about AA^ — J^N 
coupling effects in hypernuclei. It is thought that this effect might play an important role 
in neutron-rich A hypernuclei, because the total isospin becomes large. 



For such an aim, the structure of ^He, ^He, ^Li and ^Be was investigated 
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4] using an a 



cluster model, and it was pointed out that these hypernuclei were of neutron- or proton-halo 
structure. One of them, the neutron-rich A hypernucleus AHe, was observed in the (e, e'K^) 
reaction at Jlab [5], and an observed A separation energy of B\ = 5.68 ± 0.03 ± 0.25 MeV 
was reported. 

On the other hand, to produce a neutron-rich hypernucleus j^Li, the double-charge 
exchange (vr", K^) reaction on a ^°B target was performed at KEK [6|. On the basis of the 
result of this experiment, Umeya and Harada [7] calculated the structure of the Li-isotope 
A hypernuclei within a shell-model framework, and they stressed that AA^ — J^N coupling 



effects play a moderate role for the A 
Recently, a FUNUDA experiment 



Dinding energy as the number of neutrons increases. 

8|, |9|] reported an epoch-making observation of super- 
heavy hydrogen-A hypernucleus ^H as a bound state with B^ = 4.0 ±1.1 MeV. ^H is a 
neutron-rich system including four neutrons and only a single proton, which goes far toward 
the neutron drip line. Others have previously analyzed the structure of this hypernucleus. 
For example, Dalitz et al. predicted that ^H was a bound state based upon shell model 



arguments 10| • Akaishi et al. suggested that coherent AA^ — UN couphng is important for 



the A binding in ^H 11|. 



With motivation from these experimental and theoretical studies, a search (Experiment 



E-10) for ^H using the double-charge exchange (vr", A'"'") reaction at J-PARC 1^ has been 
performed, and the analysis is now in progress. Thus, it is quite timely to study the structure 
of I^H within a realistic approach. 



It should be noted that the ground state of the core nucleus H was observed as a resonant 



state with a broad width, E = l.l ± 0.3 MeV (F = 1.9 ± 0.4 MeV) [13(, with respect to 
the tnn three-body breakup threshold. We further note that the A separation energy, i?A, 
depends strongly on the spatial structure (size) of the core nucleus [lj,ll5|. Therefore, before 
calculating the binding energy of ^H, it is essential to reproduce the observed data for ^H. 
The resonant structure of ^H has been studied in the literature within the framework of 
a tnn three-body model; it was found that ^H was described well by using such a model 
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171 ]. Therefore, it is reasonable to employ a tnnK four-body model for the study of ^H 
in the present work. In order to discuss the energy and width of ^H on the basis of the 
tnn three-body model, we employ the complex scaling method (CSM), a powerful tool for 



analyzing such three-body resonance states 18|-|22[|. 



All the interactions among subunits (triton, two neutrons and A) are chosen to reproduce 
the appropriate low energy properties, such as binding energies and scattering phase shifts 
for each of the subsystems composed of two and three subunits. Using these interactions, 
we calculate the energy of ^H based on the tnnA. four-body model and discuss whether or 
not ^H exists as a bound system. 

In Sect. 2, the method used in the four-body calculation for the tnnK system is described. 
In Sect. 3, we explain the interactions employed. The calculated results and discussion are 
presented in Sect. 4. A summary is given in Sect. 5. 

II. MODEL 

A. The ^H hypernucleus 

In this work the hypernucleus ^H is considered as a triton cluster, a A particle, and 
two valence neutrons. The triton cluster is considered to be an inert core (an elementary 



particle) and to have a wave function <l>(t) with a (Os)^ shell model configuration, . 

Nine sets of Jacobian coordinates for the four-body system of ^H are illustrated in Fig.l; 
additionally, we further take into account the antisymmetrization between the two neutrons. 




FIG. 1: Jacobi coordinates for all the rearrangement channels (c = 1 ~ 9) of the tnnA four-body 
system. The two nucleons are to be antisymmetrized within each channel. 



The total Hamiltonian and the Schrodinger equation are given by 

{H-E)^jm,ttAa^) = 0, 

2 



(1) 

(2) 



i=l 



where Vtm and Vnin2 are the interactions between the triton and i-th neutron and the one 
between two neutrons, respectively. Vt^ is the triton-A interaction and Vtn^n2 is the triton- 
nn three-body force. These interactions are explained in the next section. The total wave 
function is described as a sum of amplitudes for all the rearrangement channels shown in 
Fig. [1] and LS coupling scheme: 



*Jm(aZ) - XI Yl XI XI ^nl,NL,uX,IK,ss 
c=l uLNLmX IK ss'S 



'S *(*) 
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+ X X XX ^nl,NLM\JK,ss'S ^W 
c=8 tiLNLmX IK ss'S 



Here, the operator A stands for antisymmetrization between the two valence neutrons. 
In Eq. ([3]), the isospin couphng is omitted since |H is streched in isospin space. 
Xi(^)? Xi(^)) Xi(^)) V^i''^) ^^d rii{t) are the spin and isospin functions of the neutron, 
triton and A, respectively. It is to be noted that, in Eq. ([2]), the way of spin and isospin 
coupling is the same for all the channels (c = 1 ~ 9) simply to perform numerical calculation 
easily. But, we have no problem for the present calculation, since we employ the full space 
of the spin and isopin variables for each channel. The isospin- and spin function employed 
in Eq. ([3]) are also successfully applied to the calculation of four-nucleon bound state of 



4l 



He|26|. 



Following the Gaussian Expansion Method (GEM) 23I-I25L l27l|. we take the functional 
forms of (f)nim{r), i^nlmCR) and C%{pc) as 

UM = p^e-(''/^^)V,,(p), (4) 

where the Gaussian range parameters are chosen according to geometric progressions: 

rn = ria""^ (n = 1 - ra^ax) , 

p^ = pia'"^^ (z/= 1 - Z/max) • (5) 

The angular momentum space oi l,L, X < 2 is found to be sufficient to obtain good 
convergence of the calculated results. Application of the GEM to the three-, four- and five- 
body calculations of single- and double-A hypernuclei have been extensively performed and 



are reviewed in Refs. 27H30|. 



The eigenenergies E in Eq. (1) and the coefficients C in Eq. (3) are determined by 
diagonalizing the Hamiltonian in Eq. (2) with the use of the four-body basis functions 
introduced above. If the calculated lowest state is obtained below the lowest-lying \B. + n + n 
threshold, it is identified as the the bound ground state of |H. 



When the lowest state is not a bound state but a resonant state, we calculate the energy 

n 

and width of the resonance employing the stabilization (real-scaling) method [31'] that is 
useful in calculating narrow resonances and is tractable even in the case of four-body resonant 
states. (Note that we take the CSM when calculating the broad three-body resonance of ^H 
as explained below). In the stabilization method, we first diagonalize the Hamiltonian and 
obtain the eigenenergies in the same manner as in the bound-state calculation but changing 
(scaling) the Gaussian range parameters (5), for example, as r„ — )■ arn{a ~ 0.5 to 2). 



The position and width of a resonance can be estimated, using Eq. (4) of Ref. 3l|, from 



the degree of the stabilization of the eigenenergy versus the scaling of the Gaussian range 



parameters (5). A good example of such a calculation is shown in Ref. 32|] for the study of 
the pentaquark resonances under the scattering boundary condition. 

B. The ^H nucleus 

We assume that the ^H nucleus is represented by the tnn three-body system, which is 
a part of the entire tnnA four-body system. The sets of the Jacobi coordinates describing 
the ^H nucleus are given by the channels c = 2 and c = 5 in Fig. 1 but omitting the A 
particle and the coordinate p. The wave function of ^H is written with the use of Eq. (^ 
but employing the channels c = 2 and c = 5 only and omitting the amplitudes for the A 
particle. The Hamiltonian of the tnn subsystem, say Htnn, is given by Eq. ([2]) without the 
terms for the A particle. 

As mentioned before, the ground state of ^H is a resonant state which was observed at 
Ej: = 1.7 ± 0.3 MeV with F = 1.9 ± 0.4 MeV with respect to the tnn three-body breakup 
threshold. In order to perform a study of such a broad, low-lying three-body resonance, we 



employ the CSM 18|-|22|. The CSM and its application to the nuclear physics problems are 



extensively reviewed in Ref. 33|] and references therein. Very recently, various types of 3a 



resonances in ^^C were studied using the CSM 34 1. 

Using the CSM, one can directly obtain the energy Er and the decay width F of the tnn 
three-body resonance by solving the eigenvalue problem for the complex scaled Schrodinger 
equation with a scaling angle 9: 

[Htnn{e)-E{e)]^tnn{0)=O, (6) 



where the boundary condition of the three-body outgoing wave is automatically satisfied for 
the resonance, giving E = Er — iV /2 that is, in principle, independent of 9. The complex 
scaled Hamiltonian Htnn{d) is obtained by making the radial coordinate transformation 

re -^ r^ e\ i?c -^ Re e'^ (c = 2 and 5 in Fig. 1) (7) 

in the Hamiltonian Hmn of the tnn system which is introduced in Sect. Illb. 

A great advantage of the CSM is that the resonance states are described with an L^- 
integrable wave function. Therefore, the resonance wave function '^tnn can be expanded in 
terms of the basis functions such as Eqs. (4) and (5). 

III. INTERACTIONS 



For the A^A^ interaction Vnn-, we take the AV8' potential 35|. For the AA^ interaction. 



Vaat, we employ an effective single- channel interaction simulating the basic features of the 



Nijmegen model NSC97f 36|], where the AA^-SA^ coupling effects are renormalized into AA^- 



AA^ components: The potential parameters of the AA^ interaction are listed in Table 1(a) 
in Ref. 1^ with the three-range Gaussian potentials which are chosen to reproduce the AA^ 
scattering phase shifts calculated from the NSC97f. Then, their second-range strengths in 
the even states of spin-independent and spin-spin terms are adjusted so that the calculated 
energies of the O"'"-!"'' doublet states of ^H in the NNNA four-body calculation reproduce 
the observed energies of the states. More importantly, this A — A^ interaction was applied 
to ^He within an an A three-body model, resulting in A separation energy B^ = 4.21 MeV 
which reproduces the observed B\ = 4.18 ± 0.10 MeV. This same AA^ interaction and the 
present A^A^ interaction were already applied in the calculation of the binding energy of ^He 
within the framework of an aANN four-body model, resulting in 5a = 5.48 MeV (see the 
result using 'even' in Fig.2 of Ref. [^.) which is consistent within the error bar with the 
recent data of ^He at Jlab|5|]. We thus consider that the employed AA^ interactions in this 
work are reasonable. 

The interaction VtA is obtained by folding the AA^ G-matrix interaction derived from the 
Nijmegen model F(NF) |32] into the density of the triton cluster, its strength being adjusted 
so as to reproduce the experimental value of 5a(aH) within the tA cluster model. Also, the 
spin dependence of the tA interaction is such that it reproduces the 0+-1+ doublet splitting 



in ^H. The potential parameters are listed in Table 1(b) in RefJ^ |. 

As for the Vtn, we employ a potential proposed in Ref. 38|; it is of Gaussian shape 
with dependence on the angular momentum and spin of the tn system. The d- and /-wave 
components were taken to be the same as the s and p waves, respectively (this was noted 



afterwards in Ref. 



39|). Additionally, we assume the same dependence for the partial waves 



with / > 4 (though their effect must be negligible), and describe the resultant Vtn potential 
in the parity-dependent way: 

^ r 1 + -Pr 

[ {H,even,i(^) + Vss,even,i{f) §„ " Stj 



Vt 



tn 



Yl [ {H,even,i(r) + Vss,even,i{'^) ^n " St} 
+ {Vo,oddA'^) + VssMdA'^) Sn " St} 



2 
1-Pr 



+ ^so,.(r)€-(S„ + St 



1- Pr. 



-/J'ir'- 



(8) 



where Pr is the space exchange (Majorana) operator. Sn =(Tn/2 and St=crt/2, namely, spin 
operators for the neutron and triton cluster, respectively. The potential parameters of Vtn 
are listed in Table I. 

TABLE I: Parameters of the triton-neutron potential Vtn defined by Eq. ([8]). The parameters are 
the same as in Ref. [38( except that the components for the partial waves / > 4 are omitted there. 
Size parameters are in fm~^ and strengths are in MeV. 



1 



^j 


2.1 


3.725 


3.015 


<'0,even,i 


205 


- 


- 


*^ss,even,i 


60 


- 


- 


Vb,odd,i 


- 


-2.2 


-13.95 @ 


'ss,odd,i 


- 


8.8 


-18.6 


Vso,^ 


- 


- 


4.67 



However, since the observed energy of the ground state of ^H cannot be reproduced by 
using the two-body tn potential only, we introduce an effective three-body force, Vtnn in 
Eq. ([2]), whose definition is the same as in Ref. |l7j : 



V 



tnn 



V^sfee-^/'H/^^i,)^^ 



(9) 



where the hyperradius pu is defined by 

(m„ + m„ + mt) p^ = ninrln + "^t^^ + ^tr'^n 

with the obvious notation. We take a range parameter of 636 = 2.6 fm and tune the strength 
parameter V^b so as to reproduce the energy and width of ^H within the error bar of the 
experimental data. The details are discussed in the next section. 

IV. RESULTS AND DISCUSSION 

Before studying the structure of the hypernucleus ^H, it is essentially important to 
succeed in explaining the structure of the core nucleus ^H to which a A particle is added. 
We solved the CSM equation (6) and obtained the ^H ground state as a resonant pole with 
J = 1/2+. Taking the strength of the three-body potential as V^i, = —45 MeV, we have 
a resonance pole at E,. = 1.60 MeV and F = 2.44 MeV close to the central value of the 
observed energy with the error {Ej- = 1.7 ± 0.3 MeV and F = 1.9 ± 0.4 MeV), which is 
illustrated in Fig. 2(a) for ^H (upper). 

For ^H, due to the AA^ spin-spin interaction, we have 0+-1+ spin-doublet states and the 
O"*" state is the ground state. We discuss how the energy of the 0+ ground state of ^H changes 
with respect to the position of the ^H resonance as the strength V^h of the tnn three-body 
force is increased. 

As shown in Fig. 2(a), when the energy of ^H is close to the central value of the observed 
data, the ground state of ^H is obtained, with V^b = —45 MeV, at the total energy E = —0.87 
MeV (F = 0.23 MeV), which is located by 1.13 MeV above the {R + n + n threshold as 
a resonance. In this case, we have a A separation energy B\ = 2.47 MeV with respect to 
^H + A. 

It is expected that the 1"*" state of ^H is located above the O"*" state by about 1 MeV. This 
is because the two valence neutrons have spin in ^H, and the AA^ spin-spin interaction 
between the two neutrons and the A particle cancels. The energy splitting of the O^-l"*" 
doublet states in ^H is considered to correspond to that in ^H. Therefore, in the case of (a) 
in Fig. |2l the 1^ state is expected to lie above the tnnA four-body breakup threshold. At 
present, it is difficult to precisely calculate the energy and decay width of the 1^ resonance 
using the stabilization method. 



1/2^ (E=i-66,r=2.4 4)" ----- "|Exp. E=1.7±0.3 

"(a) "(b)" 

(E=1.17,r=0.91) 

OMeV 

t+n+n 



t+n+n+A 
OMeV 



(a)(E=-0.87,r=0.23) 



^H+n+n 

-2.00 
(b) -2.07 



FIG. 2: Calculated energies and decay widths of ^H (upper) and ^H (lower) in the case that the 
tnn three-body force (9) has a strength (a) V^i, = —45 MeV and (b) V^i, = —65 MeV. 

As the energy of ^H is lowered by increasing the strength of the three-body force, the 
energy of ^H becomes closer to the lowest ^H + n + n threshold. However, even if we adjust 
the calculated energy of ^H to the lower edge of the error band on the observed energy of 
^H, we cannot obtain any bound state for ^H. 

When the resonance energy of ^H reaches E^ = 1.17 MeV with V^b = —65 MeV, we 
obtain a very weakly bound state of ^H at E = —2.07 MeV as shown in Fig. Mjo). In this 
case we have B\ = 3.24 MeV with respect to ^H + A, which is consistent with the FINUDA 



data 5a = 4.0 ± 1.1 MeV 



0.9 



within the experimental error. 
In order to demonstrate how Bj>^{\}i) depends on the energy and decay width of the core 
nucleus ^H, we consider one more case. If we have a tnn three-body force with Vsb = —73 
MeV, then we obtain a lower-lying resonance for ^H at E,. = 0.73 MeV with F = 0.50 MeV 
and a bound state of ^H at E = -3.02 MeV that is 1.01 MeV below the XH + n + n threshold. 
In this case, we have B\ = 3.75 MeV, which is much larger than B\ = 2.47 MeV in the case 
of (a) in Fig. |2l Here, we note that the energy and width of the above artificially lowered 
^H resonance is similar to those of the ground-state resonance of the ^He nucleus [E,. = 0.80 
MeV with F = 0.65 MeV), but B^ of ^H (3.75 MeV) is smaller than Ba of ^He (4.18 ± 0.10 



MeV by experiment and 4.21 MeV by our calculation). This fact is reasonable. Because, 
two valence nucleons in ^H are occupied in the Op orbit. On the other hand, in ^He one is 
in the Op orbit and the other is in Os orbit. Then, this results in the weaker AA^ attraction 
(smaller Ba) in ^H than in ^He . 

As mentioned before, even if we adjust the tnn three-body force so as to match the 
lower limit of the error of the observed energy of ^H, we could have no bound state for ^H. 
However, by identifying \R ground state through the observation of its two-body weak-decay 
7r^ meson, the FINUDA experiment J8|, |9[ was able to claim that ^H is a bound system. 
Then, in our tnnA four-body model with no explicit AA^-SA^ coupling effect, the ^H resonant 
state should exist above the tnn three-body breakup threshold by less than 1.17 MeV. 

Recently, the E-10 experiment to search ^H has been performed at J- PARC and the 
analysis is now in progress. If the experiment confirms the existence of this hypernucleus as 
a bound state, it is desirable to measure the binding energy of the core nucleus ^H with a 
precision of 100 keV. 

Let us compare our results for the A separation energy B\ with those given in the 
literature. The shell model analysis by Dalitz 10|] reported B\ = 4.2 MeV. A recent shell 
model calculation by Millener cited in Ref. |9| gave Ba = 4.28 ± 0.10 MeV by using the 
three binding energies of iHe {Ba = 3.12 ± 0.02 MeV), lHe(SA = 5.36 ± 0.09 MeV), and 
\B.{Ba = 2.04 ± 0.04 MeV). The results for Ba of ^ H from those shell model calculations 
are similar to the observed data of ^H. Akaishi et al ll|], however, obtained Ba = 5.8 MeV, 
which is by ~ 1.5 MeV larger than the above two shell- model results. Each of these three 
results for Ba are much larger than the Ba from the present calculation. The difference comes 
from whether or not the core nucleus ^H is taken into account explicitly in the calculations 
as a broad three-body resonant state. 

Here, it should be reiterated that AA^ — SA^ coupling is not explicitly treated in our 
model. Since the total isospin of ^H is 3/2, it is likely that the coupling plays an important 
role, especially, working as an effective ANN three-body force in the binding energies of ^H. 
So far, some authors have investigated the role of the AA^ — J^N coupling in the binding 



energies of ^H, ^H and y^^He 40|-|47|. For example, in Ref. 46|, four-body calculations for 



^H and \B.e were performed taking account of the AN-TiN coupling explicitly, and it was 
found that the effective ANN three-body force contributes about 0.6 MeV to the binding 
energy of the O"*" ground state in the case of the NSC97f YN interaction. As far as the 



n 

neutron-rich hypernucleus ^H is concerned, Akaishi et al. [11] argued that KN-TjN couphng 
gave a large contribution to its binding energy. On the other hand, a recent shell model 



calculation by Millener cited in Ref. [48| found that the AA^ — SA^ contribution is small in 
^H. Thus, the determination of the size of the effect in the binding energy of ^H is still 
unsettled. To answer the question, it is necessary to perform a coupled-channel calculation 
taking into account the tnnK and t(^He)A^A^S channels. This is one of our future subjects. 
The new data from the search for ^H in Experiment E-10 at J- PARC will provide us with 
useful information about AA^-SA^ coupling. 

Finally, it is interesting to examine the spatial structure of the ground state of ^H 
when the state becomes bound as in Fig. 2(b). Using our wave function of this state, we 
calculated the single-particle density distributions of the constituent particles and plotted 
them in Fig. |3l the dotted curve denotes the density of the Os nucleon in the triton, the 
dashed curve is for the A particle and the solid line for one of the two valence neutrons. 
The A particle is bound to the triton cluster mostly in the Osa orbit that is much shallower 
than the Os orbit of the nucleon in the triton, whereas the valence nucleons are very loosely 
coupled to the triton. Thus, we see, as should be anticipated, that there are three layers 
in the matter distribution of the hypernucleus ^.H, namely, the nucleons in the well bound 
triton core, the loosely bound A skin in the ^H subsystem, and the two-neutron halo that 
is unbound in ^H. 

V. SUMMARY 

We have studied the structure of the neutron rich hypernucleus ^H within the framework 
of a tnnK four-body model. For this study of ^H, it is essentially important to reproduce the 
property of the ground state of the core nucleus ^H which is a low-lying, broad three-body 
resonance at E^ = 1.7±0.3 MeV (F = 1.9±0.4 MeV) with respect to the t + n + n threshold. 
We note that the A separation energy, Sa, depends strongly on the spatial size of the core 
nucleus [M, [l5|. We thus utilize a tnn three-body model for the ^H nucleus and treat both 
^H and ^.H consistently. 

In the present tnnA model, all the two-body interactions among subunits (triton, two 
neutrons and A) are chosen to reproduce the low energy properties, such as binding energies 
and scattering phase shifts of each of the subsystems composed of two and three subunits. 




FIG. 3: Single-particle densities of the Os nucleons comprising the triton (dotted curve), the A 
particle (dashed curve) and one of the valence neutrons (solid curve) in the ground state of the 
hypernucleus ^H when the state becomes a bound state shown in Fig. 2(b). The radius r is 
measured from the cm. of the triton. 



The NN interaction is given by the AV8' potential 35| . We employ a AA^ effective potential 
that simulates AA^ scattering phase shifts of the NSC97f interaction 36|] and is slightly tuned 
to reproduce the observed energies of the spin-doublet 0^-1"'" states of ^H. The AA^ — UN 
coupling is renormalized into the AA^ potential. But, we note that there remains an effective 
ANN three-body force that is not renormalizeable into the effective AA^ potential. The 
employed AA^ interaction reproduces the observed binding energies of ^He and ^Hej^] within 
the framework of aAA^ and aANN three- and four-body models, respectively. 



38 



39[. We found that the observed resonance 



The tn potential is taken from Ref. 
energy of ^H cannot be reproduced with the two-body A^A^ and tn interactions only, and 
therefore we introduced a phenomenological, attractive tnn three-body force. When the tnn 
force is tuned to reproduce the central value of the observed ^H resonance energy, we obtain 
Ej. = 1.60 MeV (r = 2.44 MeV) and, at the same time, we have the ^H ground state as a 
resonance at E = —0.87 MeV (F = 0.23 MeV) with respect to the t + n + n + A threshold; 
it is located 1.13 MeV above the lowest ^H + n + n threshold (Fig. 2(a)). In this case, we 
have i?A = 2.47 MeV. Even if the tnn three-body force were adjusted to reproduce the lower 
band {E,. = 1.4 MeV) based on the error of the observed ^H resonance energy, we could not 



obtain any bound state for ^H below the ^H + n + n threshold. 

We were able to generate a shallow bound state of ^H at E = —2.07 MeV (by 0.07 
MeV below the ^H + n -\- n threshold) if the tnn three-body force is tuned to have the ^H 
resonance at E = 1.17 MeV (F = 0.91 MeV) which is, however, 0.23 MeV below the quoted 
lower band of the error on the observed energy (Fig. 2(b)). In this case, we have B\ = 3.24 
MeV, which is consistent with the FINUDA data {B\ = 4.0 ± 1.1 MeV) within the error. 

In order to study the structure of ^H more fully, we anticipate having a more precise 
value for the resonance energy of ^H (100 keV accuracy). It should be noted that we did 
not explicitly include AA^ — UN coupling in the present work. If, in the analysis of the E-10 
experiment at J-PARC, now in progress, ^H is confirmed to exist as a bound state, we could 
obtain useful information about AA^ — SA^ coupling. To study the effect of such coupling, it 
is essential to perform a coupled tnnA + t(^He)A^A^S four-body calculation. This calculation 
will be one of our future endeavors. 
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